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Abstract.  This  paper  evaluates  the  performance  of  vari¬ 
ous  linear  ozone  photochemistry  parametrizations  using  the 
stratosphere-troposphere  data  assimilation  system  of  the  Met 
Office.  A  set  of  experiments  were  run  for  the  period 
23  September  2003  to  5  November  2003  using  the  Cari¬ 
olle  (vl.O  and  v2.1),  LINOZ  and  Chem2D-OPP  (vO.l  and 
v2.1)  parametrizations.  All  operational  meteorological  ob¬ 
servations  were  assimilated,  together  with  ozone  retrievals 
from  the  Michelson  Interferometer  for  Passive  Atmospheric 
Sounding  (MIPAS).  Experiments  were  validated  against  in¬ 
dependent  data  from  the  Halogen  Occultation  Experiment 
(HALOE)  and  ozonesondes.  Additionally,  a  simple  offline 
method  for  comparing  the  parametrizations  is  introduced. 

It  is  shown  that  in  the  upper  stratosphere  and  mesosphere, 
outside  the  polar  night,  ozone  analyses  are  controlled  by  the 
photochemistry  parametrizations  and  not  by  the  assimilated 
observations.  The  most  important  factor  in  getting  good  re¬ 
sults  at  these  levels  is  to  pay  attention  to  the  ozone  and  tem¬ 
perature  climatologies  in  the  parametrizations.  There  should 
be  no  discrepancies  between  the  climatologies  and  the  assim¬ 
ilated  observations  or  the  model,  but  there  is  also  a  compet¬ 
ing  demand  that  the  climatologies  be  objectively  accurate  in 
themselves.  Conversely,  in  the  lower  stratosphere  outside  re¬ 
gions  of  heterogeneous  ozone  depletion,  the  ozone  analyses 
are  dominated  by  observational  increments  and  the  photo¬ 
chemistry  parametrizations  have  little  influence. 

We  investigate  a  number  of  known  problems  in  LINOZ 
and  Cariolle  vl.O  in  more  detail  than  previously,  and  we  find 
discrepancies  in  Cariolle  v2.1  and  Chem2D-OPP  v2.1,  which 
are  demonstrated  to  have  been  removed  in  the  latest  avail- 
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able  versions  (v2.8  and  v2.6  respectively).  In  general,  how¬ 
ever,  all  the  parametrizations  work  well  through  much  of  the 
stratosphere,  helped  by  the  presence  of  good  quality  assimi¬ 
lated  MIPAS  observations. 


1  Introduction 

Atmospheric  ozone  is  of  interest  not  just  because  of  the  prob¬ 
lem  of  ozone  depletion  (e.g.,  WMO,  2003)  but  also  for  its  ba¬ 

sic  role  in  controlling  the  temperature  structure  of  the  atmo¬ 
sphere  through  the  absorption  of  solar  and  long-wave  radia¬ 
tion.  The  distribution  of  atmospheric  ozone  results  from  the 
interaction  of  transport,  chemistry,  and  radiation  processes. 
A  full  description  of  the  photochemistry  of  ozone  would 
be  extremely  complicated,  involving  hundreds  of  chemical 
species  and  reactions,  many  of  which  are  interlinked.  A  de¬ 
tailed  approximation  to  this  can  be  embodied  in  a  chemistry- 
transport  model  (CTM)  (e.g.  Chipperfield,  1999;  Rozanov 
et  ah,  1999;  Errera  and  Fonteyn,  2001;  Josse  et  al.,  2004), 
which  is  driven  by  off-line  meteorological  analyses.  How¬ 
ever,  in  applications  such  as  climate  modelling  and  numer¬ 
ical  weather  prediction  (NWP)  or  data  assimilation  (DA),  it 
is  often  desirable  to  use  a  faster,  simpler  representation  of 
ozone  photochemistry. 

Cariolle  and  Deque  (1986)  developed  a  parametrization 
of  ozone  photochemistry  based  on  a  linearisation  of  the 
ozone  tendency  around  an  equilibrium  state,  using  param¬ 
eters  derived  from  a  CTM  with  more  detailed  chemistry. 
The  parametrization  depends  only  on  temperature  and  ozone 
amount;  hence  no  other  chemically  active  species  need 
be  modelled.  As  well  as  saving  computer  resources,  this 
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prevents  a  mis-specified  or  poorly-known  chemical  species 
from  causing  a  bias  in  the  modelled  ozone  amount.  Along 
with  the  original  scheme  of  Cariolle  and  Deque  (1986,  vl.O), 
an  updated  version  has  recently  become  available  (Cariolle 
and  Teyssedre,  2007,  v2),  and  coefficients  have  been  inde¬ 
pendently  developed  by  McLinden  et  al.  (2000,  the  LINOZ 
scheme),  McCormack  et  al.  (2004,  Chem2D-OPP  vO.l)  and 
McCormack  et  al.  (2006,  Chem2D-OPP  v2).  Table  1  sum¬ 
marises  these  developments.  Linearised  ozone  photochem¬ 
istry  schemes  are  routinely  used  in  data  assimilation  (e.g., 
Riishpjgaard,  2000;  Eskes  et  al.,  2003;  Dethof  and  Holm, 
2004;  Geer  et  al.,  2006a,b)  and  in  multi-year  climate  sim¬ 
ulations  (e.g.,  Hadjinicolaou  et  al.,  1997).  Recently,  Taylor 
and  Bourqui  (2005)  developed  a  fast  ozone  photochemistry 
scheme  of  intermediate  complexity,  and  this  may  be  very 
useful  in  the  future.  Historically,  one  of  the  simplest  ap¬ 
proaches  has  been  to  make  ozone  loss  rates  proportional  to 
the  ozone  amount  and  to  keep  production  rates  fixed.  This 
was  tried  in  the  GEOS  ozone  data  assimilation  system,  but 
major  limitations  were  found  (Riishpjgaard,  2000;  McCor¬ 
mack  et  al.,  2006).  This  paper  examines  only  the  much-used 
linear  approach  of  Cariolle  and  Deque  (1986). 

A  number  of  studies  have  indicated  problems  in  the  earlier 
generation  of  linear  photochemistry  schemes  and  the  ozone 
distributions  that  they  produce.  LINOZ  is  unsuitable  for  use 
in  the  upper  stratosphere  (McCormack  et  al.,  2004;  Geer 
et  al.,  2006a).  Version  1.0  of  Cariolle  and  Deque  (1986) 
has  a  strong  sensitivity  to  the  overlying  ozone  amount  (Geer 
et  al.,  2006b).  Even  small  ozone  differences  can  be  impor¬ 
tant  in  general  circulation  models  (GCMs):  ~10%  variations 
in  ozone  amounts  can  result  in  changes  in  modelled  tempera¬ 
tures  of  several  Kelvin  (Cariolle  and  Morcrette,  2006).  Now 
that  updated  photochemistry  schemes  are  available,  it  is  use¬ 
ful  to  evaluate  and  understand  their  differences.  The  data 
assimilation  (DA)  framework  presents  an  opportunity  to  do 
this. 

We  examine  the  performance  of  a  number  of  differ¬ 
ent  ozone  photochemistry  schemes  in  the  stratosphere- 
troposphere  DA  system  of  the  Met  Office.  3-D-Variational 
data  assimilation  (3D-Var)  is  used  to  assimilate  all  opera¬ 
tional  dynamical  observations,  plus  ozone  retrievals  from  the 
Michelson  Interferometer  for  Passive  Atmospheric  Sounding 
(MIPAS)  on  the  Envisat  satellite.  The  analyses  are  validated 
against  independent  ozone  data  from  sondes  and  the  Halogen 
Occultation  Experiment  (HALOE). 

Using  a  data  assimilation  system  to  evaluate  a  number  of 
different  parametrizations  is  a  relatively  novel  approach.  It 
is  more  typical  to  do  such  evaluations  with  multi-year  free 
model  runs.  However,  free-  running  models  can  evolve  to 
a  state  that  may  be  different  from  the  real  atmosphere;  the 
regular  insertion  of  observational  data  in  a  DA  system  acts 
to  prevent  this.  With  DA,  model  forecasts  or  analyses  can 
be  compared  directly  to  observations  in  their  synoptic  con¬ 
text.  However,  there  are  limitations.  A  normally  useful  prop¬ 
erty  of  DA  is  that  slow-growing  model  errors  are  swiftly  cor¬ 


rected,  but  here,  this  means  that  such  errors  cannot  easily  be 
detected.  Also,  the  NWP  system  used  here  is  computation¬ 
ally  expensive  and  experiments  could  be  run  only  for  short 
periods.  A  better  choice  in  future  for  such  experiments  would 
be  a  CTM-based  data  assimilation  system,  computationally 
much  cheaper.  This  paper  does  not  eliminate  the  need  to 
examine  linear  ozone  parametrizations  over  long  time  peri¬ 
ods  using  model  free-runs.  Instead  it  provides  a  validation  in 
the  type  of  system  in  which  such  parametrizations  are  most 
likely  to  be  used,  and  it  uncovers  a  number  of  issues  which 
are  specific  to  DA.  However,  it  does  also  identify  a  number  of 
fast-growing  errors  in  some  schemes,  and  these  are  important 
to  both  DA  and  free-running  models. 


2  Ozone  chemistry  schemes 

2.1  Photochemistry  parametrization 


The  rate  of  change  of  ozone  due  to  photochemistry  can  be 
written  (Cariolle  and  Deque,  1986;  McLinden  et  al.,  2000; 
McCormack  et  al.,  2006)  as  a  first  order  Taylor  expansion 
about  the  ozone  production  rate  ( P )  minus  loss  rate  (L),  at 
an  equilibrium  state; 


C  =  (P  —  L)  o  + 


d(P-L) 

dX 


(X  -  Xo)  + 

0 


d(P  ~  L) 
dT 


(T  -  To)  + 
o 


HP  -  L) 
9$ 


(<b  -  O0). 
o 


(1) 


This  parametrization  has  three  variables:  the  local  ozone 
mixing  ratio,  /  (in  this  paper  as  mass  mixing  ratio  in  units  of 
kg  kg-1);  the  temperature,  T ,  in  K,  and  the  column  of  ozone 
overlying  the  level  under  consideration,  <t>,  in  kgm  ,  where; 


l 

®=1  /g  f 

TOA 


X  dp. 


(2) 


and  the  integral  runs  over  all  pressure  levels  from  the  top  of 
the  atmosphere  (TOA)  down  to  level  l,  where  /  is  the  level 
under  consideration.  All  other  items  in  Eq.  (1)  are  coeffi¬ 
cients  valid  at  the  equilibrium  state  (denoted  with  the  sub¬ 
script  0)  which  are  either  climatological  values  or  have  been 
pre-calculated  with  a  detailed  photochemical  model.  They 
are  given  as  functions  of  latitude,  model  level  and  month; 
hence  there  is  no  diurnal  or  longitudinal  variation.  These  in¬ 
clude  the  equilibrium  ozone  production  rate  minus  loss  rate, 
(P—L)o  and  its  partial  derivatives,  and  the  climatological 
values  of  ozone,  /o>  temperature,  7q,  and  the  overlying  col¬ 
umn  of  ozone,  <J>o,  which  is  calculated  from  vertical  profiles 
of  xo  using  Eq.  (2). 

The  second  term  in  the  expansion  accounts  for  variations 
in  the  local  ozone  amount,  the  third  for  temperature  and  the 
last  term  for  the  influence  of  non-local  ozone  on  the  amount 
of  solar  radiation  reaching  the  level  in  question,  and  here  we 
will  call  it  the  radiation  term. 
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Table  1.  Linear  ozone  photochemistry  coefficients. 


Name 

Notes 

Reference 

Cariolle  vl.O 

Cariolle  and  Deque  (1986) 

Cariolle  vl.2 

As  vl.O  but  with  heterogeneous 
chem.  term 

Cariolle  and  Deque  (1986);  Dethof  and  Holm 
(2004) 

Cariolle  v2 

Cariolle  and  Teyssedre  (2007) 

LINOZ 

McLinden  et  al.  (2000) 

Chem2D-OPP  vO.l 

(P—L)  o  and  ozone  terms  only 

McCormack  et  al.  (2004) 

Chem2D-OPP  v2 

All  four  terms 

McCormack  et  al.  (2006) 

This  photochemistry  parametrization  is  primarily  intended 
for  use  in  the  stratosphere.  It  causes  no  diurnal  variation, 
but  it  should  be  remembered  that  there  is  in  reality  a  diurnal 
cycle  of  ozone  which  starts  to  become  significant  at  levels 
above  0.5  hPa,  in  the  mesosphere.  Also,  the  lack  of  longitu¬ 
dinal  variation  is  not  appropriate  in  the  troposphere,  where 
the  surface  sources  of  ozone  precursors  show  strong  tempo¬ 
ral  and  spatial  heterogeneity. 

In  the  DARC  implementation,  the  radiation  term  has  been 
modified  by  substituting  the  following  into  Eq.  (1): 

l 

T>  -  =  \/g  f  (x  -  xo  )dp.  (3) 

TOA 


We  have  used  Eqs.  (2)  and  (3)  to  eliminate  the  overlying  col¬ 
umn  climatology,  <J>0i  from  the  parametrization.  Though  it 
may  seem  trivial,  this  modification  means  that  <J>  and  <J>o 
are  calculated  implicitly  using  consistent  methods.  If  they 
were  not  calculated  consistently,  the  radiation  term  could 
produce  a  forcing  even  if  x  —  Xo  throughout  the  overlying 
column.  With  this  modification,  <t>o  no  longer  needs  to  be 
recalculated  from  xo  each  time  the  model’s  vertical  resolu¬ 
tion  changes.  Neglecting  to  do  this,  and  instead  interpolating 
<$>o  to  a  new  set  of  vertical  levels,  caused  a  ~40%  model  bias 
in  the  lower  mesosphere  in  the  DARC/Met  Office  assimila¬ 
tion  experiments  examined  in  Geer  et  al.  (2006a). 

We  can  explore  the  effect  of  the  parametrization  by  consid¬ 
ering  the  ozone  budget  in  a  hypothetical  model,  which  will 
include  the  rates  of  change  (or  tendency)  of  ozone  due  to 
modelled  ozone  transport,  A ,  and  chemistry,  C  (from  Eq.  1): 

3x 

-7T  =  A  +  C.  (4) 

at 

If  the  temperature,  overhead  column  ozone  and  transport  (T , 
<l>  and  A  respectively)  are  constant  then  we  can  define  a 
steady  state  mixing  ratio: 


Xss  =  Xo+ 


(P-L)  o  + 


B(P-L) 
3  T 


(T  -  T0)  + 
o 


d(P  -  L) 
3$ 


(<D  -  O0)  +  A 
o 


r, 


and  the  photochemical  relaxation  time, 


(5) 


r  =  -1/ 


(d(P-L)  \ 

V  3X  o/  ' 


(6) 


This  allows  us  to  solve  Eq.  (4)  analytically,  showing  that  with 
steady-state  transport  and  constant  temperature  and  radia¬ 
tion  terms,  the  linearised  ozone  photochemistry  scheme  will 
cause  modelled  ozone  to  follow  an  exponential  relaxation  to 
steady  state,  Xss >  with  a  time  constant  r: 

Xt+A‘  =  x'  +  (Xss  -  x')d  -  e“Af/T),  (7) 


Ozone  amounts  x‘  and  X,+Al  apply  at  times  t  and  t+At  re¬ 
spectively.  Equations  (5),  (6)  and  (7)  are  similar  to  those 
derived  in  McLinden  et  al.  (2000)  and  McCormack  et  al. 
(2006).  However,  the  steady  state  (Eq.  5)  is  different  here, 
because  the  effects  of  transport  have  been  included. 

Figure  1  shows  the  October  value  of  r,  the  photochemi¬ 
cal  relaxation  time,  from  the  Cariolle  vl.O  parametrization. 
In  the  upper  stratosphere,  outside  the  polar  night,  r  is  less 
than  a  day  and  the  ozone  fields  are  essentially  in  photo¬ 
chemical  equilibrium.  If  perturbed  away  from  this  equilib¬ 
rium,  ozone  fields  will  quickly  relax  back  to  it.  Lower  in 
the  stratosphere,  photochemical  relaxation  times  are  much 
longer,  with  r>100  days  at  100  hPa.  Here,  modelled  trans¬ 
port  has  an  important  control  over  the  ozone  field. 

When  temperature  and  overhead  column  ozone  are  in  their 
climatological  state,  i.e.  T—Tq  and  the  steady  state 

mixing  ratio  can  be  simplified: 


Xss  =  Xo  +  [(T5  -  L) o  +  A]r. 


(8) 


Equation  (8)  encapsulates  the  climatological  balance  be¬ 
tween  photochemistry  and  transport.  If  the  steady  state 
ozone  amount  generated  by  the  photochemistry  scheme  is 
to  be  equal  to  climatological  ozone,  xo-  Eq.  (8)  requires 
either  that  modelled  ozone  transport  must  balance  the  net 
equilibrium  production  or  loss  due  to  photochemistry  (i.e. 


www.atmos-chem-phys.net/7/939/2007/ 


Atmos.  Chem.  Phys.,  7,  939-959,  2007 


942 


Fig.  1.  October  values  of  r,  the  photochemical  relaxation  time,  in 
days,  from  the  Cariolle  vl.O  parametrization. 

(P—L)o+A—0),  or  that  r->0.  Hence  in  this  example,  ozone 
will  always  be  close  to  climatology  in  regions  where  the  pho¬ 
tochemical  relaxation  time  is  short.  In  the  mid  and  lower 
stratosphere  and  troposphere,  relaxation  times  are  relatively 
long  (Fig.  1).  Here,  the  steady-state  example  reflects  real-life 
behaviour:  In  the  lower  stratosphere,  for  example,  net  pho¬ 
tochemical  production  in  the  tropics  is  approximately  bal¬ 
anced  by  the  upward  and  poleward  transport  of  ozone  in  the 
Brewer-Dobson  circulation  (e.g..  Plumb  et  al.,  2002);  there 
is  a  net  photochemical  loss  at  higher  latitudes. 

2.2  Practical  considerations 

In  the  case  of  the  radiation  term,  we  have  already  seen  that  an 
incorrectly  implemented  linear  photochemistry  scheme  can 
generate  a  spurious  photochemical  tendency  (Eq.  3  and  fol¬ 
lowing  discussion).  The  same  principle  applies  to  the  other 
terms  in  the  photochemistry  scheme. 

Equations  (7)  and  (8)  show  that  if  modelled  ozone  trans¬ 
port  A  does  not  balance  the  equilibrium  photochemical  ten¬ 
dency  (P— L)q,  modelled  ozone  will  relax  toward  a  steady 
state  ( Xss )  that  is  different  from  climatology  ('/<))■  If  we  de¬ 
note  the  imbalance  as  an  ozone  tendency  e  then  that  steady 
state  would  be: 

Xss  =  Xo  +  ex.  (9) 

If  modelled  ozone  transport  were  different  from  climatology 
because  of  natural  variability,  of  course  this  would  be  a  de¬ 
sirable  result.  However,  particularly  in  data  assimilation  sys¬ 
tems,  stratospheric  constituent  transport  can  be  erroneously 
fast  (e.g.,  Schoeberl  et  al.,  2003).  Equation  (9)  shows  that 
errors  in  the  ozone  tendency,  or  in  transport,  have  most  effect 
on  modelled  ozone  where  the  photochemical  relaxation  time 
r  is  relatively  long,  i.e.  in  the  mid  and  lower  stratosphere. 
Geer  et  al.  (2006b)  showed  how  an  excessively  fast  modelled 
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Brewer-Dobson  circulation  caused  biases  in  the  DARC/Met 
Office  ozone  analyses  in  the  mid  stratosphere. 

In  the  troposphere,  it  is  important  to  account  for  ozone 
transport  by  convection  and  boundary  layer  processes.  In  the 
CTM  used  to  generate  the  photochemistry  coefficients,  these 
processes  may  be  parametrized  through  the  use  of  vertical 
diffusion.  In  the  GCM  in  which  the  ozone  photochemistry 
parametrization  is  to  be  used,  convection  and  boundary  layer 
processes  will  be  more  explicitly  resolved,  though  still  heav¬ 
ily  parametrized.  If  modelled  transport  A,  differs  between 
the  CTM  and  the  GCM,  the  equilibrium  photochemical  ten¬ 
dency  (P—L) o  will  be  consistent  with  the  CTM,  but  not  the 
GCM.  Again,  this  could  cause  ozone  amounts  to  be  biased  in 
the  GCM. 

In  data  assimilation  applications  it  is  necessary  that  model 
and  observations  should  be  unbiased  with  respect  to  one  an¬ 
other.  If  the  model  were  relaxing  to  steady  state  (Eq.  7), 
would  this  state  be  unbiased  with  respect  to  the  climatology 
of  assimilated  ozone  observations?  It  is  usual  to  find  small 
biases  between  different  instruments  and  different  climatolo¬ 
gies.  Hence,  it  is  typical  in  data  assimilation  applications 
(e.g.,  Eskes  et  al.,  2003)  to  use  a  new  ozone  climatology  (/o) 
in  place  of  the  one  supplied  with  the  photochemistry  scheme. 
It  is  hoped  that  the  new  climatology  will  be  less  biased  with 
respect  to  observations.  Often,  as  in  this  study,  the  Fortuin 
and  Kelder  (1998)  climatology  is  chosen.  Later  we  suggest 
that  it  may  actually  be  necessary  to  use  a  climatology  based 
on  the  observations  being  assimilated. 

Similar  considerations  apply  to  the  temperature  climatol¬ 
ogy.  If  climatological  mean  temperatures  in  the  model  were 
different  from  those  used  in  the  photochemistry  scheme.  To, 
the  temperature  term  (see  Eq.  1)  would  on  average  pro¬ 
duce  a  net  forcing  in  ozone,  with  a  consequent  effect  on  the 
steady  state  (Eq.  5).  Even  if  climatological  transport  and 
photochemistry  were  balanced  ((P— L)o+A—0),  and  over¬ 
head  column  ozone  was  in  its  equilibrium  state  (<J>=<t>o)-  the 
model  would  relax  to  a  steady  state  that  would  be  different 
from  climatology.  If  the  erroneous  forcing  were  denoted  e 
then  Eq.  (9)  would  again  apply. 

Figure  2  compares  the  monthly  mean  analysed  temper¬ 
ature  for  October  2003  to  To  from  three  of  the  chemistry 
schemes  used  here.  The  climatologies  and  the  DARC/Met 
Office  analyses  show  differences  of  up  to  20  K.  Between  6 
and  30  hPa  at  80°  S  there  is  a  warm  bulge  in  Met  Office  tem¬ 
peratures,  compared  to  climatology,  and  this  is  likely  due 
to  the  strong  minor  warming  that  took  place  during  Octo¬ 
ber  2003  (Lahoz  et  al.,  2006,  see  also  Fig.  13).  Randel 
et  al.  (2004)  found  that  the  CIRA86  climatology,  as  used  in 
Chem2D-OPP  v2.1,  has  a  5-10  K  warm  bias  through  much 
of  the  stratosphere.  This  is  consistent  with  the  positive  bias 
seen  in  CIRA86  temperatures  in  Fig.  2,  compared  to  the  other 
climatologies,  particularly  around  the  stratopause. 

The  results  later  in  this  paper  illustrate  the  problems  that 
can  be  caused  by  discrepancies  between  the  parametriza- 
tion’s  ozone  (/o)  and  temperature  (To)  climatologies  and  the 
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Fig.  2.  October  mean  temperature  at  (a)  80°  S,  (b)  40°  S,  (c)  the  Equator,  from  the  October  2003  DARC/Met  Office  analyses  (solid),  the 
Cariolle  vl.O  (dot-dashed)  and  v2.1  (dashed)  coefficients,  and  CIRA  1986  climatology  used  in  Chem2D  OPP  v2.1  (dotted). 


modelled  or  observed  equivalents.  However,  a  more  careful 
choice  of  xo  and  Tq  will  likely  improve  analysed  ozone  dis¬ 
tributions.  In  contrast,  there  is  no  easy  solution  to  the  prob¬ 
lem  of  erroneous  or  mis-matched  ozone  transport;  this  par¬ 
ticularly  affects  modelled  lower  stratosphere  ozone  amounts. 
Data  assimilation  can  correct  the  ozone  distribution  here,  but 
not  the  underlying  model  errors. 

2.3  Heterogeneous  chemistry 


Heterogeneous  ozone  chemistry  must  be  modelled  in  order 
to  describe  ozone  depletion  in  the  spring  polar  vortex  (e.g. 
WMO,  2003).  In  all  experiments,  heterogeneous  chemistry 
is  parametrized  by  a  cold  tracer  scheme  similar  to  that  in¬ 
troduced  by  Hadjinicolaou  et  al.  (1997).  The  focus  here  is 
on  photochemistry,  and  the  heterogeneous  chemistry  is  kept 
fixed. 

The  cold  tracer  parametrization  used  here  (P.  Braesicke, 
personal  communication)  imposes  an  additional  forcing  on 
the  ozone  field  to  simulate  the  heterogeneous  depletion  of 
ozone; 


dx 

dt 


het 


(10) 


Here,  i]  is  a  time  constant  for  heterogeneous  ozone  depletion, 
set  to  10  days,  and  b  is  the  cold  tracer,  an  advected  3-D  field. 
The  forcing  is  applied  only  in  direct  sunlight  (solar  zenith  an¬ 
gle  <90°).  The  cold  tracer  represents  the  degree  of  chlorine 
activation;  its  chemical  evolution  is  described  by: 


db 


dt 


chem 


i  i 

—  (1  -b) - b. 

1 IP  t]L 


(ID 


The  first  term  on  the  right  hand  side  is  active  only  below 
a  temperature  threshold  representative  of  polar  stratospheric 
cloud  (PSC)  formation,  Tp$c,  and  describes  the  rapid  ac¬ 
tivation  of  heterogeneous  chemistry  with  a  time  constant 
rjp—  4h.  Ip  sc  varies  between  191.5  K  at  26.5  km  and  203  K 
at  11.5  km.  The  term  is  switched  off  outside  this  altitude 
range,  and  equatorward  of  55°.  The  second  term  is  active 
only  in  direct  sunlight  and  describes  the  slow  loss  of  chlorine 
activation  with  a  time  constant  of  in—\Q  days  in  the  southern 
hemisphere  (SH)  and  r]p—  5  days  in  the  northern  hemisphere 
(NH).  Equation  (11)  implies  that  the  cold  tracer  b  is  bounded 
by  0  (representing  no  chlorine  activation)  and  1  (full  activa¬ 
tion).  In  practice,  the  cold  tracer  shows  nearly  full  chlorine 
activation  in  the  SH  polar  vortex  in  the  lower  stratosphere 
in  our  experiments  for  October  2003,  with  ozone  depletion 
taking  place  only  in  sunlit  areas. 

A  similar  version  of  this  scheme  has  been  used  succes- 
fully  in  the  KNMI  (Royal  Netherlands  Meterological  Insti¬ 
tute)  ozone  data  assimilation  system  (Eskes  et  al.,  2003)  and 
was  able  to  make  a  good  simulation  of  ozone  depletion  even 
without  data  assimilation  (Siegmund  et  al.,  2005).  The  As¬ 
similation  of  Envisat  Data  (ASSET)  intercomparison  (Geer 
et  al.,  2006a)  examined  ozone  analyses  based  on  models  with 
many  different  treatments  of  heterogeneous  chemistry,  in¬ 
cluding  DARC/Met  Office  ozone  analyses  made  with  a  sys¬ 
tem  similar  to  that  used  here.  The  DARC/Met  Office  anal¬ 
yses  performed  adequately  well  in  the  ozone-hole,  though 
amounts  were  not  depleted  to  the  near-zero  values  observed 
by  sondes.  This  was  due  not  to  deficiencies  in  the  cold  tracer 
scheme,  but  instead  to  erroneous  ozone  production  in  the 
radiation  term  of  the  Cariolle  vl.O  photochemistry  scheme. 
Later  we  show  that  ozone  depletion  is  well  represented  in  the 
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analyses  when  the  cold  tracer  scheme  is  used  with  any  of  the 
other  photochemistry  parametrizations. 

Cariolle  and  Teyssedre  (2007)  have  recently  described  two 
parametrizations  of  heterogeneous  chemistry  that  may  be 
used  with  v2  of  their  photochemistry  scheme.  The  first  one 
is  a  simple  extra  term  in  Eq.  (1),  active  in  sunlight  at  tem¬ 
peratures  less  than  195  K,  describing  heterogeneous  ozone 
depletion  as  proportional  to  the  ozone  amount.  The  second 
is  a  variant  of  the  cold  tracer  technique.  We  do  not  test  either 
of  those  parametrizations  here,  but  Cariolle  and  Teyssedre 
(2007)  show  that  each  parametrization  is  capable  of  repre¬ 
senting  polar  ozone  depletion  in  multi-annual  model  runs, 
with  the  cold  tracer  providing  a  longer  maintenance  of  low 
ozone  in  the  vortex,  and  additional  export  of  ozone  depletion 
to  mid-latitudes.  Geer  et  al.  (2006a)  have  also  shown  that 
in  assimilation  experiments  in  the  presence  of  good  quality 
ozone  observations,  the  Cariolle  v2.1  chemistry  scheme  with 
the  simple  (non  cold-tracer)  heterogeneous  parametrization 
is  capable  of  producing  results  in  the  ozone  hole  similar  to 
much  more  detailed  treatments. 

2.4  Cariolle  and  Deque  (1986)  vl.0,  v2.1  and  v2.8 

The  Cariolle  vl.0  scheme  is  described  by  Cariolle  and  Deque 
(1986)  and  was  calculated  using  a  2-D  photochemical  model 
with  an  upper  boundary  at  1  hPa,  and  by  extrapolation  or 
from  1-D  model  results  above  that  level. 

The  subsequent  v2  scheme  is  described  by  Cariolle  and 
Teyssedre  (2007).  It  has  been  derived  using  the  same  2-D 
photochemical  model  as  for  vl.0,  though  there  have  been 
a  number  of  changes.  Gas-phase  chemical  rates  have  been 
upgraded  using  the  JPL  (2003)  recommendations,  and  to¬ 
tal  chlorine  is  set  to  year  2000  amounts.  In  contrast,  vl.0 
is  based  on  knowledge  of  chemical  rates  and  constituent 
amounts  from  the  early  1980s.  Another  change  in  v2  is  that 
the  temperature  distribution  and  the  residual  meridional  cir¬ 
culation,  used  for  minor  tracer  transport,  are  derived  from  a 
10  year  simulation  of  the  Arpege-Climat  general  circulation 
model.  This  was  found  to  be  as  important  to  the  resulting 
coefficients  as  the  change  in  chemical  rates.  The  vertical  and 
horizontal  resolutions  of  the  2-D  model  have  also  been  in¬ 
creased  to  match  the  Arpege-Climat  discretisation.  The  2-D 
model  has  45  vertical  levels  extending  up  to  0.1  hPa,  64  lati¬ 
tudes,  and  accounts  for  the  photochemistry  of  63  species,  29 
of  which  are  transported. 

The  partial  derivatives  appearing  in  Eq.  (1)  are  obtained 
by  perturbing  the  2-D  model  fields  by  ±10%  for  the  ozone 
mixing  ratio  (more  precisely  the  odd  oxygen  family)  and  the 
ozone  column,  and  by  ±10  K  for  the  temperature.  For  each 
perturbed  case  the  non-transported  short  lived  species  are  re¬ 
evaluated  at  steady  state  and  the  resulting  ozone  production 
and  loss  rates  are  used  in  the  partial  derivatives  calculations. 
This  is  done  for  every  month  and  a  set  of  7  zonal  mean  coef¬ 
ficients  are  obtained. 
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In  our  assimilation  experiments  we  tested  v2.1  of  the  Car¬ 
iolle  scheme.  Close  to  publication,  we  were  able  to  include 
the  latest  version,  v2.8,  but  only  in  the  offline  comparisons 
presented  in  Sect.  2.7.  There  we  see  that  v2.8  has  improved 
accuracy  in  the  upper  stratosphere  and  mesosphere,  and  in 
the  lower  troposphere.  Three  configurations  of  the  Cariolle 
parameters  are  tested  with  assimilation  experiments:  (a)  us¬ 
ing  the  supplied  ozone  climatology  in  the  v2.1  scheme,  (b) 
instead  using  Fortuin  and  Kelder  (1998)  climatology  in  the 
v2.1  scheme,  and  (c)  using  Fortuin  and  Kelder  (1998)  clima¬ 
tology  in  the  vl.0  scheme. 

Note  that  neither  of  the  Cariolle  v2  heterogeneous  ozone 
depletion  parametrizations  are  tested  here;  the  cold  tracer 
scheme  described  in  Sect.  2.3  was  used  in  all  experiments. 

2.5  LINOZ 

The  LINOZ  scheme  is  described  by  McLinden  et  al.  (2000). 
The  coefficients  in  Eq.  (1)  were  calculated  for  12  months  be¬ 
tween  85°  S  and  85°  N  and  and  at  25  altitudes  between  10  to 
58  km  using  a  photochemical  box  model  (Prather  and  laffe, 
1990;  Prather,  1992).  This  model  includes  109  kinetic  re¬ 
actions,  36  photolysis  reactions  and  43  species,  and  reaction 
coefficients  and  absorption  cross-sections  are  adopted  from 
DeMore  et  al.  (1997).  Species  concentrations  are  character¬ 
istic  of  the  1990s  (Avallone  and  Prather,  1997).  For  each 
monthly  calculation,  the  box  model  is  integrated  for  30  days 
with  the  diurnal  cycle  fixed  at  mid-month.  The  ozone  ten¬ 
dency  and  partial  derivatives  (see  Eq.  1)  are  diurnally  aver¬ 
aged. 

Originally,  the  ozone  and  column  ozone  climatologies 
from  McPeters  (1993)  and  the  temperature  climatology  from 
Nagatani  and  Rosenfield  (1993)  were  used.  Here  instead 
we  substitute  the  ozone  climatology  of  Fortuin  and  Kelder 
(1998),  consistent  with  the  way  LINOZ  was  used  in  ozone 
data  assimilation  experiments  at  KNMI  (Eskes  et  al.,  2003). 

LINOZ  coefficients  are  not  available  below  10  km.  Geer 
et  al.  (2006a)  found  that  in  the  troposphere,  a  relaxation 
to  ozone  climatology  produced  smaller  biases  compared  to 
ozonesonde  than  did  the  full  linear  chemistry  parametriza¬ 
tions.  Hence  for  LINOZ  in  these  tests,  below  10  km,  a  re¬ 
laxation  to  Fortuin  and  Kelder  (1998)  climatology  was  im¬ 
plemented,  using  the  photochemical  relaxation  times  of  Car¬ 
iolle  vl.0. 

2.6  Chem2D-OPP  vO.l,  v2.1  and  v2.6 

Chem2D-OPP  is  described  by  McCormack  et  al.  (2004), 
for  vO.l,  and  by  McCormack  et  al.  (2006)  for  v2;  see 
http://uap-www.nrl.navy.mil/dynamics/html/chem2dopp/ 
chem2d_opp.html  for  updates.  Photochemistry  coefficients 
are  computed  with  the  NRL-Chem2D  middle  atmosphere 
photochemical-transport  model  (Siskind  et  al.,  2003).  The 
Chem2D  model  domain  extends  from  pole  to  pole  and 
from  the  surface  up  to  p—2  x  10-5  hPa  (~122km  altitude). 
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with  47  vertical  levels.  The  middle  atmospheric  radiative 
heating,  photochemistry,  and  transport  are  fully  coupled. 
Chem2D  photochemistry  accounts  for  54  chemical  species 
and  uses  updated  JPL  (2003)  reaction  rates.  Halogen  species 
amounts  are  constant  in  the  troposphere  and  are  taken  from 
WMO  (2003);  stratospheric  amounts  are  model-determined 
after  a  20  year  spin-up  with  the  troposphere  as  a  boundary 
condition.  Diurnally  averaged  photolysis  rates  are  computed 
by  averaging  hourly  values,  and  diurnally  averaged  reaction 
rate  coefficients  are  derived  from  pre -computed  night-day 
ratios  of  relevant  species.  Solar  irradiance  for  calculating 
both  photolysis  and  UV  radiative  heating  is  specified  as 
a  function  of  wavelength  from  120-800  nm,  including 
Lyman-a  and  both  the  Schumann-Runge  continuum  and 
Schumann-Runge  bands  for  O2  photolysis. 

Chem2D-OPP  vO.l  includes  the  first  two  terms  from  the 
right  hand  side  of  Eq.  (1)  and  does  not  account  for  tempera¬ 
ture  or  radiation  effects  (McCormack  et  al.,  2004).  Diurnally 
averaged  values  of  the  Chem2D  net  ozone  tendency  (P—L) 0 
are  computed  for  the  15th  day  of  each  month.  The  partial 
derivative  with  respect  to  mixing  ratio,  31  1  |o  in  Eq.  (1), 

is  calculated  as  the  negative  inverse  of  the  ozone  photochem¬ 
ical  relaxation  time  r.  The  latter  is  determined  from  the  sum 
of  individual  loss  rates  involving  reactions  with  NOx,  Clx, 
and  HOx.  The  results  are  tabulated  as  functions  of  latitude, 
pressure,  and  month  and  then  interpolated  onto  a  10°  latitude 
grid  at  standard  pressure  levels  from  1000-0.001  hPa. 

Chem2D-OPP  v2.1  includes  all  four  of  the  terms  in 
Eq.  (1).  The  v2.1  ( P—L)0  and  3(  |o  coefficients  are 
computed  as  in  Chem2D-OPP  vO.  1 .  To  evaluate  the  temper¬ 
ature  and  column  ozone  coefficients,  |o  and  —  Pd^L)  |o, 

respectively,  the  Chem2D  model  computes  (P—L) 0  for  a 
given  altitude,  latitude,  and  time  of  year,  then  immediately 
repeats  the  calculation  using  identical  model  constituent 
fields  and  a  perturbation  in  either  temperature  or  overlying 
ozone  column  amount.  In  the  temperature  case,  perturba¬ 
tions  (AT)  between  ±20  K  are  imposed  and  the  entire  chem¬ 
ical  system  is  then  solved  with  an  iterative  Newton-Raphson 
technique  until  the  (P—L)  values  converge  to  a  new  equilib¬ 
rium  state.  Similarly,  the  coefficient  31  '  jo  is  evaluated  by 

introducing  ozone  column  perturbations  AO  between  ±50% 
to  the  Chem2D  UV  transmission  functions  used  to  compute 
the  O2  and  O3  photolysis  rates.  For  a  more  detailed  descrip¬ 
tion  of  this  method  see  McCormack  et  al.  (2006). 

The  Chem2D  model  uses  fixed  heating  rates  at  the  surface 
as  a  model  boundary  condition  and  so  the  radiative  heating 
is  not  coupled  to  model  ozone  in  the  lowermost  model  lev¬ 
els.  For  implementation  in  the  DARC  system,  the  v2.1  radi¬ 
ation  term  was  turned  off  below  500  hPa  as  it  was  suspected 
it  would  not  work  well  in  the  lower  troposphere.  Surface 
values  of  the  v2.1  photochemical  relaxation  time  were  er¬ 
roneously  negative  due  to  an  error  in  the  vertical  interpola¬ 
tion  scheme,  and  this  would  have  caused  a  runaway  growth 
in  ozone  amounts.  To  prevent  this  happening,  r  was  reset 
to  ~+2  days  at  the  surface.  The  experiments  in  this  paper 
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Fig.  3.  Key  to  the  parametrizations  shown  in  Figs.  4  to  11. 


also  indicate  a  number  of  other  discrepancies  in  v2.1.  As  for 
the  Cariolle  scheme,  we  were  able  at  a  late  stage  to  examine 
the  most  recent  Chem2D-OPP  coefficients  (v2.6),  but  again 
only  in  the  offline  comparisons.  These  comparisons  show 
that  most  discrepancies  have  now  been  removed.  Updated 
values  for  3( Pd~L)  |o,  which  now  include  Brx  effects,  provide 
shorter  relaxation  times  in  the  lower  stratosphere,  bringing 
them  closer  to  the  Cariolle  schemes. 


2.7  Offline  comparison  of  photochemistry  schemes 


We  can  examine  the  relative  strengths  of  the  terms  in  the  dif¬ 
ferent  schemes  by  testing  their  sensitivity  to  a  representative 
perturbation  in  ozone  or  temperature.  For  ozone,  we  used  a 
perturbation  based  on  the  climatological  ozone  standard  de¬ 
viations  of  Fortuin  and  Kelder  (1998),  given  as  a  function  of 
month,  pressure  and  latitude.  For  the  overlying  ozone  col¬ 
umn,  we  calculated  the  partial  column  integral  of  these  stan¬ 
dard  deviations  using  Eq.  (2).  For  temperature,  we  assumed 
a  uniform  perturbation  of  5  K,  though  around  the  wintertime 
polar  vortex,  stratospheric  temperatures  can  vary  by  much 
larger  amounts. 

For  each  month,  latitude  and  pressure  level,  the  change  in 
the  net  photochemical  ozone  tendency  caused  by  the  pertur¬ 
bation  was  normalised  by  the  climatological  ozone  amount, 
Xo-  For  example,  the  sensitivity  to  changes  in  local  ozone 
was  calculated  as: 


AC  = 


<rx  d(P  -  L) 
XO  3x 


0 


(12) 


where  ax  is  the  climatological  ozone  standard  deviation.  We 
then  converted  AC  to  units  of  %  per  day.  The  net  clima¬ 
tological  ozone  tendency,  (P—L)q,  was  normalised  by  the 
climatological  ozone  amount  so  that  it  could  be  examined 
in  the  same  units.  The  resulting  sensitivities  are  shown  for 
the  months  of  January,  April,  July  and  October  at  the  50hPa 
level  (Fig.  4)  and  the  5  hPa  level  (Fig.  5). 

It  is  clear  from  Figs.  4  and  5  that  though  the  sensitivity 
of  the  coefficients  can  vary  both  zonally  and  from  month  to 
month,  the  major  differences  between  the  parametrizations 
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Fig.  4.  Rate  of  change  of  ozone  (%  per  day)  produced  by  the  photochemistry  schemes  at  the  50hPa  level,  from  the  (P—L) q  term  (a-d)  and 
in  response  to  typical  perturbations  of  ozone  (e-h),  temperature  (i— 1),  and  overlying  column  ozone  (m-p).  Figures  are  shown  for  January  (a, 
e,  i,  m),  April  (b,  f,  j,  n),  July  (c,  g,  k,  o)  and  October  (d,  h,  1,  p).  See  colour  key  in  Fig.  3.  Note  only  one  line  is  shown  for  Cariolle  v2.1, 
because  these  values  are  independent  of  the  climatology  coefficients  (in  later  figures  we  need  to  distinguish  which  ozone  climatology  was 
used).  Note  also  that  Chem2D-OPP  vO.l  has  no  temperature  or  radiation  term,  so  it  shows  zero  sensitivity  in  those  figures. 


Atmos.  Chem.  Phys.,  7,  939-959,  2007 


www.atmos-chem-phys.net/7/939/2007/ 


Radiation  term:  Temperature  term:  Ozone  term:  (P-L)0term: 

effect  on  ozone  /%  per  day  effect  on  ozone  /%  per  day  effect  on  ozone  /%  per  day  effect  on  ozone  /%  per  day 


A.  J.  Geer  et  al.:  Linear  ozone  photochemistry  parametrizations 


947 


a  January 


April 


July 


d  October 


m 


Latitude 


Latitude 


Latitude 


Latitude 


Fig.  5.  As  for  Fig.  4,  except  at  the  5  hPa  level.  Note  that  in  panels  (m-p)  the  Chem2D-OPP  v2. 1  and  v2.6  curves  are  superposed. 
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Fig.  6.  Global  mean  absolute  values  of  the  rate  of  change  of  ozone  (%  per  day)  produced  by  the  October  coefficients  of  the  photochemistry 
schemes  from  the  (P—L) o  term  (a)  and  in  response  to  typical  perturbations  of  ozone  (b),  temperature  (c),  and  overlying  column  ozone  (d). 
See  colour  key  in  Fig.  3.  Note  that  zero  rates  of  change  cannot  be  shown  on  this  figure  because  of  the  logarithmic  scale;  for  example  there  is 
no  line  associated  with  the  Chem2D-OPP  v2.1  radiation  term  at  500  hPa  and  below,  where  this  term  is  set  to  zero. 


are  relatively  consistent  over  all  latitudes  and  months.  To 
summarise  Figs.  4  and  5  still  further,  and  to  extend  this  anal¬ 
ysis  to  all  vertical  levels,  we  have  calculated  global  means 
(equal  weight  by  latitude)  from  the  absolute  value  of  AC, 
for  each  pressure  level  and  month.  Figure  6  shows  these 
global  mean  sensitivities  for  the  month  of  October;  other 
months  are  very  similar  (not  shown).  Because  the  figure 
shows  global  mean  absolute  sensitivities,  there  are  no  neg¬ 
ative  values.  In  general,  the  schemes  have  least  influence  on 
the  ozone  amount  around  the  tropopause,  and  larger  influ¬ 
ence  in  the  troposphere,  upper  stratosphere  and  mesosphere, 
where  ozone  photochemistry  is  faster.  Of  course,  the  relax¬ 
ation  times  in  Fig.  1  show  a  similar  picture  for  the  ozone 
term,  but  we  are  now  able  to  compare  between  all  terms  in 
the  parametrization. 

Figures  4a-d,  5a-d  and  6a  show  that  outside  the  tropical 
lower  stratosphere  the  (P—L) o  term  is  small  compared  to  the 
response  of  the  other  terms  to  representative  perturbations. 
With  the  exception  of  LINOZ,  and  outside  the  lower  tropo¬ 
sphere,  the  magnitude  of  the  (P—L) o  term  is  similar  in  all 
schemes.  At  50hPa  (Figs.  4a-d),  representative  the  lower 
stratosphere,  there  is  a  net  ozone  production  in  the  trop¬ 
ics  and  net  destruction  in  the  higher  latitudes,  which  should 
largely  be  balanced  by  ozone  transport  in  the  Brewer-Dobson 
circulation  (see  Sect.  2.1).  In  the  tropics,  net  ozone  pro¬ 


duction  dominates  over  the  other  terms  of  the  parametriza¬ 
tion.  At  5  hPa  (Figs.  5a-d,  excepting  LINOZ)  the  effect  of 
the  (P—L) o  term  hovers  around  zero  net  production.  Here, 
and  even  more  strongly  at  higher  levels  in  the  stratosphere, 
the  other  terms  of  the  parametrization  dominate  over  the 
(P—L) o  term. 

It  is  thought  there  may  still  be  small  discrepancies  in 
our  knowledge  of  photochemistry  in  the  upper  stratosphere, 
given  the  observed  ~10%  ozone  difference  between  mod¬ 
els  and  observations  (e.g.  Natarajan  et  al.,  2002).  However, 
even  if  these  discrepancies  were  to  result  in  a  small  error  in 
(P—L) o  in  the  parametrizations  tested  here,  its  importance 
would  be  small:  First  because  the  sensitivities  of  the  other 
terms  in  the  parametrization  are  much  larger  (Fig.  6),  and 
second  because  any  error  would  be  limited  in  its  effect  by 
the  short  relaxation  times  at  these  levels  (Eq.  9). 

Above  lOhPa,  the  LINOZ  (P—L) o  term  is,  erroneously, 
orders  of  magnitude  larger  than  in  the  other  schemes,  and 
hence  this  has  an  effect  despite  the  short  relaxation  times  at 
these  levels.  Above  lOhPa,  in  LINOZ,  this  term  causes  a 
strong  net  loss  of  ozone,  so  large  negative  biases  are  gen¬ 
erated  above  lOhPa  (McCormack  et  al.,  2004;  Geer  et  al., 
2006a). 

Another  likely  error  is  that  the  (P—L) o  terms  in  both 
Cariolle  vl.O  and  v2.1  are  excessively  large  in  the  lower 
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troposphere,  being  of  order  10%  per  day  (Fig.  6a).  No  ac¬ 
curacy  is  claimed  for  these  schemes  in  the  troposphere,  but 
they  should  at  least  be  well  behaved.  Later  we  see  that  such 
excessive  net  production  results  in  relatively  large  ozone  bi¬ 
ases  in  the  lower  troposphere. 

Figures  4e-h,  5e-h  and  6b  show  the  sensitivity  of  the  dif¬ 
ferent  parametrizations  to  ozone  variations.  To  aid  under¬ 
standing,  Fig.  7  shows  the  corresponding  relaxation  times 
(Eq.  6)  at  the  equator.  Note  that  there  is  no  direct  equivalence 
between  the  graphs:  Figs.  4  to  6  depend  on  latitude-varying 
factors  as  well  as  the  reciprocal  of  the  relaxation  time  (see 
Eq.  12).  The  apparent  discrepancy,  for  example  at  5hPa  at 
the  equator  in  October,  between  relaxation  times  of  2  days 
(Fig.  7)  and  sensitivities  of  less  than  2%  per  day  (Fig.  5h),  is 
explained  by  the  climatological  standard  deviation,  which  is 
only  3%  here. 

Both  the  Cariolle  vl.O  and  v2.1  schemes  are  up  to  a  fac¬ 
tor  of  10  more  sensitive  than  LINOZ  and  Chem2D-OPP 
vO.l  and  v2.1  to  ozone  variations  in  the  lower  stratosphere 
and  at  the  tropopause.  In  terms  of  photochemical  relaxation 
time,  this  corresponds  to  r~100  days  at  the  tropopause  com¬ 
pared  to  r  approaching  1000  days  in  LINOZ  and  Chem2D- 
OPP.  This  difference  was  noted  by  McCormack  et  al.  (2004), 
though  it  is  now  thought  that  the  rapid  changes  in  ozone  they 
found  in  Cariolle  vl.O  hindcast  experiments  at  high  north¬ 
ern  latitudes  in  the  mid-stratosphere  (around  lOhPa,  their 
Fig.  11),  which  were  attributed  to  the  shorter  relaxation  time, 
were  in  fact  most  likely  due  to  the  very  strong  radiation  term 
in  Cariolle  vl.O. 

Version  2.6  of  the  Chem2D-OPP  coefficients  includes  cat¬ 
alytic  cycles  involving  bromine  compounds  (Brx),  producing 
somewhat  shorter  values  of  r  in  the  lower  stratosphere.  We 
see  that  these  values  remain  slightly  longer  than  correspond¬ 
ing  ones  in  the  Cariolle  vl.O  scheme.  It  must  be  noted  that 
the  Cariolle  schemes  and  Chem2D-OPP  use  different  meth¬ 
ods  to  compute  the  relaxation  time.  In  Cariolle’s  schemes  it 
is  computed  after  allowing  for  readjustment  of  the  concentra¬ 
tions  of  short  lived  species  in  response  to  the  ozone  pertur¬ 
bation,  whereas  Chem2D-OPP  takes  an  instantaneous  value. 

In  the  middle  and  upper  stratosphere  where  the  ozone  pro¬ 
duction  is  dominated  by  the  photodissociation  of  O2  the  two 
methods  should  converge.  However,  the  two  approaches  may 
differ  in  the  lower  stratosphere  where,  for  example,  readjust¬ 
ments  in  the  amount  of  NOx  species  will  have  a  significant 
effect  on  ozone  production.  The  validity  of  each  approach 
would  depend  on  the  timescales  of  the  perturbations:  for  fast 
perturbations  the  instantaneous  approach  should  be  valid;  for 
perturbations  with  timescales  longer  than  a  day,  the  readjust¬ 
ment  of  minor  species  should  be  taken  into  account. 

In  the  upper  stratosphere  and  mesosphere,  Cariolle  v2.1 
has  relatively  low  sensitivity  to  typical  ozone  variability, 
compared  to  the  other  schemes.  Equivalently,  this  means  that 
photochemical  relaxation  times  are  substantially  longer  in 
Cariolle  v2.1.  This  low  sensitivity  (extending  also  to  the  tem¬ 
perature  and  radiation  terms)  resulted  from  problems  with 
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Fig.  7.  October  values  of  r,  the  photochemical  relaxation  time,  at 
the  equator,  for  each  of  the  different  parametrizations  shown  in  the 
key  in  Fig.  3. 

making  semi-implicit  calculations  in  regions  with  fast  chem¬ 
istry.  The  problem  was  fixed  by  using  an  explicit  calculation, 
and  we  see  that  the  resulting  Cariolle  v2.8  is  very  close  to  the 
majority  of  the  other  schemes. 

Chem2D-OPP  v2.1  shows  problems  near  the  poles  in 
April  and  October  (Figs.  5f  and  h).  This  was  due  to  cor¬ 
rupted  outputs  from  the  Chem2D  model,  and  has  later  been 
corrected:  Chem2D-OPP  v2.6  now  agrees  with  the  majority 
of  the  other  schemes.  However,  this  error  is  a  main  factor  in 
causing  problems  in  the  south  polar  upper  stratosphere  in  our 
assimilation  experiments  with  Chem2D-OPP  v2. 1 . 

Figures  4i— 1,  5i-l  and  6c  show  that  there  are  few  large 
discrepancies  between  the  temperature  terms  of  the  differ¬ 
ent  schemes,  though  differences  up  to  an  order  of  magnitude 
can  be  found  in  some  areas. 

Figure  4m-p,  5m-p  and  6d  show  that  in  the  radiation  term, 
the  main  outlier  is  the  Cariolle  vl.O,  which  is  excessively 
strong  compared  to  other  schemes.  It  is  not  clear  why  the 
Cariolle  vl  .0  term  should  be  so  strong,  since  similar  methods 
were  used  to  create  vl.O  and  v2.  Since  the  vl.O  coefficients 
were  distributed  many  years  ago  it  is  not  now  possible  to 
re-examine  this  in  detail.  The  strong  Cariolle  vl.O  radiation 
term  was  responsible  for  the  problem  of  excessive  creation 
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of  ozone  in  the  ozone  hole  described  in  Geer  et  al.  (2006b) 
and  the  rapid  ozone  changes  see  in  hindcast  experiments  at 
10  hPa  by  McCormack  et  al.  (2004,  their  Fig.  11).  The  other 
schemes  also  show  differences,  typically  of  less  than  an  order 
of  magnitude,  but  we  see  that  the  influence  of  the  radiation 
term  is  typically  smaller  than  that  of  the  ozone  or  temperature 
terms,  so  these  differences  are  not  particularly  important. 

In  summary,  the  LINOZ  (P—L) o  term  and  the  Cariolle 
vl.O  radiation  term  have  sensitivities  up  to  two  orders  of 
magnitude  different  from  the  equivalent  terms  in  the  other 
parametrizations.  These  erroneous  coefficients  lead  to  prob¬ 
lems  that  have  been  seen  in  a  number  of  studies  (McCor¬ 
mack  et  al.,  2004;  Geer  et  al.,  2006a,b).  In  the  early  ver¬ 
sions  of  the  most  recent  parametrizations,  Cariolle  v2.1  and 
Chem2D-OPP  v2.1,  there  were  discrepancies  of  up  to  1  order 
of  magnitude  in  some  areas,  compared  to  the  other  schemes. 
These  discrepancies  were  in  part  due  to  erroneous  calcula¬ 
tions,  but  also  in  the  lower  stratosphere,  relaxation  times  be¬ 
came  shorter  in  Chem2D-OPP,  and  closer  to  Cariolle  v2,  af¬ 
ter  Brx  chemistry  was  included.  The  latest  schemes  (Cariolle 
v2.8  and  Chem2D-OPP  v2.6)  show  only  relatively  minor  re¬ 
maining  differences,  likely  resulting  from  the  varied  ways  in 
which  the  coefficients  have  been  derived. 

3  Method 

3.1  Assimilation  system 

The  Met  Office  NWP  system  has  recently  been  extended 
to  allow  the  assimilation  of  ozone  (Jackson  and  Saunders, 
2002;  Jackson,  2004)  but  ozone  is  not  assimilated  opera¬ 
tionally.  Here,  MIPAS  v4.61  ozone  is  assimilated  in  re¬ 
analysis  mode,  alongside  all  operational  dynamical  observa¬ 
tions,  using  a  stratosphere/troposphere  version  of  the  oper¬ 
ational  NWP  system.  The  system  is  that  described  in  Geer 
et  al.  (2006a),  except  that  MIPAS  temperatures  are  no  longer 
assimilated,  since  it  was  found  that  their  assimilation  could 
degrade  analysed  stratopause  temperatures.  The  GCM  has  a 
horizontal  resolution  of  3.75°  longitude  by  2.5°  latitude  and 
50  levels  in  the  vertical,  from  the  surface  to  ~0. 1  hPa.  It  uses 
a  new  dynamical  core  (Davies  et  al.,  2005)  which  includes  a 
semi-Lagrangian  transport  scheme.  The  ozone  tracer  is  sub¬ 
ject  to  convective  and  boundary  layer  transport.  There  is  no 
feedback  between  ozone  and  radiation:  heating  rate  calcula¬ 
tions  are  done  using  an  ozone  climatology.  Data  assimila¬ 
tion  uses  3D-Var  (Lorenc  et  al.,  2000).  Ozone  is  assimilated 
univariately,  but  3D-Var  does  not  infer  dynamical  informa¬ 
tion,  so  the  only  effect  of  ozone  on  the  dynamical  analysis  is 
through  its  influence  on  temperature  radiance  assimilation. 

3.2  MIPAS 

MIPAS  is  an  interferometer  for  measuring  infrared  emissions 
from  the  atmospheric  limb  (Fischer  and  Oelhaf,  1996).  MI¬ 
PAS  operational  data  are  available  between  July  2002  and 
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March  2004,  after  which  instrument  problems  meant  it  could 
only  be  used  on  an  occasional  basis.  The  operational  mea¬ 
surements  were  made  along  17  discrete  lines-of-sight  in  the 
reverse  of  the  flight  direction  of  ENVISAT,  with  tangent 
heights  between  8  km  and  68  km.  The  vertical  resolution  was 
~3  km  in  the  stratosphere  and  the  horizontal  resolution  was 
~300km  along  the  line  of  sight.  ENVISAT  follows  a  sun- 
synchronous  polar  orbit,  allowing  MIPAS  to  sample  globally, 
and  to  produce  up  to  ~  1 000  atmospheric  profiles  per  day. 
From  the  infrared  spectra,  ESA  retrieved  profiles  of  pres¬ 
sure,  temperature,  ozone,  water  vapour,  HNO3,  NO2,  CH4 
and  N2O  at  up  to  17  tangent  points  (ESA,  2004).  MIPAS 
version  4.61  data,  reprocessed  offline,  is  used  here.  When 
treated  as  a  point  retrieval,  MIPAS  ozone  has  only  small  bi¬ 
ases  when  compared  to  independent  data  except  in  the  lower 
stratosphere  (100  to  30hPa),  where  positive  biases  of  order 
10%  are  seen  (Geer  et  al.,  2006a). 

Apart  from  a  small  number  rejected  for  quality  control  rea¬ 
sons,  all  available  MIPAS  ozone  observations  were  assim¬ 
ilated  into  the  Met  Office  NWP  system  in  the  experiments 
presented  here. 

3.3  Experiments 

Table  2  summarises  the  ozone  chemistry  characteristics  of 
the  six  assimilation  experiments  performed  here,  which 
were  otherwise  identical.  Experiments  were  initialised  on 
23  September  2003  with  fields  from  the  DARC/Met  Office 
analyses  produced  for  the  ASSET  intercomparison  and  de¬ 
scribed  in  Geer  et  al.  (2006a),  and  were  run  until  5  Novem¬ 
ber  2003.  The  period  was  limited  to  six  weeks  by  the  com¬ 
putational  expense  of  running  the  assimilation  system.  The 
strength  of  these  assimilation  experiments  is  that  we  can  test 
the  parametrizations  under  conditions  of  rapid  synoptic  vari¬ 
ability.  Hence  the  period  was  chosen  for  the  rapid  variability 
linked  with  the  breakdown  of  the  SH  polar  vortex  and  the  de¬ 
velopment  of  the  NH  polar  vortex.  It  also  captures  the  time 
of  the  deepest  extent  of  the  ozone  hole.  Lahoz  et  al.  (2006) 
describe  the  2003  SH  winter  and  spring  in  more  detail. 

3.4  Validation  framework 

Analyses  are  compared  to  independent  data  from  ozonesonde 
and  HALOE,  and  also  to  the  assimilated  MIPAS  observa¬ 
tions,  using  the  methods  described  in  Geer  et  al.  (2006a). 
Analyses  are  interpolated  onto  a  set  of  fixed  pressure  lev¬ 
els  and  sampled  daily  at  00:00  Z  and  12:00  Z,  before  being 
compared  to  observations.  The  independent  observations  are 
described  briefly  below;  see  Geer  et  al.  (2006a)  for  further 
details.  It  is  worth  noting  that  MIPAS,  sonde  and  HALOE 
have  different  temporal  and  spatial  sampling.  For  example 
MIPAS  sampled  most  latitudes  daily;  HALOE  observations 
come  from  discrete  latitude  bands  (see  coverage  plots  in  Geer 
et  al.,  2006a).  Hence,  some  differences  are  to  be  expected 
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Table  2.  Summary  of  experiments. 


Name 

Ozone  climatology,  xo 

Heterogeneous  chemistry 

Cariolle  vl.O  with  F&K  climatology 

Fortum  and  Kelder  (1998) 

Cold  tracer 

Cariolle  v2.1 

Cariolle  v2.1 

Cold  tracer 

Cariolle  v2.1  with  F&K  climatology 

Fortuin  and  Kelder  (1998) 

Cold  tracer 

LINOZ 

Fortuin  and  Kelder  (1998) 

Cold  tracer 

Chem2D-OPP  vO.l 

Fortuin  and  Kelder  (1998) 

Cold  tracer 

Chem2D-OPP  v2.1 

Fortuin  and  Kelder  (1998) 

Cold  tracer 

when  comparing  to  different  datatypes,  simply  because  of 
the  varying  geographical  and  temporal  coverage. 

3.4.1  Ozonesondes 

Ozonesondes  are  used  as  independent  data  to  validate  the 
analyses.  Profiles  have  been  obtained  from  the  World 
Ozone  and  Ultraviolet  Radiation  Data  Centre  (WOUDC, 
http://www.woudc.org/).  Southern  Hemisphere  Additional 
Ozonesondes  project  (SHADOZ,  http://croc.gsfc.nasa.gov/ 
shadoz/,  Thompson  et  al.,  2003a,b)  and  the  Network  for  the 
Detection  of  Stratospheric  Change  (NDSC,  http://www.ndsc. 
ncep.noaa.gov/).  We  use  all  available  ozonesonde  ascents, 
except  for  the  Indian  sondes  as  they  have  large  errors.  No 
other  selection  criteria  were  applied.  The  data  used  come 
from  42  locations  and  were  made  using  a  variety  of  measure¬ 
ment  techniques.  Sondes  typically  make  measurements  from 
the  surface  to  around  the  10  hPa  level.  Total  error  for  the  most 
common  type  of  ozonesonde  is  estimated  to  be  within  —7% 
to  +17%  in  the  upper  troposphere,  ±5%  in  the  lower  strato¬ 
sphere  up  to  lOhPa  and  —14%  to  +6%  at  4hPa  (Komhyr 
et  al.,  1995).  Errors  are  higher  in  the  presence  of  steep  ozone 
gradients  and  where  ozone  amounts  are  low. 

3.4.2  HALOE 

HALOE  (Russell  et  al.,  1993)  is  used  as  independent  data  to 
validate  the  analyses.  HALOE  uses  solar  occultation  to  de¬ 
rive  atmospheric  constituent  profiles,  making  the  data  sparse 
in  time  and  space,  with  about  15  observations  per  day  at  each 
of  two  latitudes.  The  horizontal  resolution  is  495  km  along 
the  orbital  track  and  the  vertical  resolution  is  about  2.5  km. 
We  use  a  version  19  product,  screened  for  cloud  using  the 
algorithm  of  Hervig  and  McHugh  (1999),  and  available  from 
the  HALOE  website  (http://haloedata.larc.nasa.gov/).  Ver¬ 
sion  19  ozone  retrievals  are  nearly  identical  to  those  of  vl8, 
and  above  the  120  hPa  level  they  agree  with  ozonesonde  data 
to  within  10%  (Bhatt  et  al.,  1999).  Below  this  level,  profiles 
can  be  seriously  affected  by  the  presence  of  aerosols  and  cir¬ 
rus  clouds. 


4  Results 


Figures  8,  9  and  10  show,  respectively,  the  mean  differences 
between  analyses  and  HALOE,  sonde  and  MIPAS  observa¬ 
tions,  given  as  a  percentage  relative  to  an  ozone  climatology. 
This  ozone  climatology  is  described  in  Geer  et  al.  (2006a) 
and  combines  those  of  Fortuin  and  Kelder  (1998)  and  Logan 
(1999).  Statistics  are  calculated  for  the  period  27  September 
2003  to  5  November  2003. 

In  the  troposphere,  upper  stratosphere  and  mesosphere, 
and  polar  vortex,  there  are  differences  between  experiments. 
Here,  biases  between  analyses  and  independent  data  can  in 
general  be  attributed  to  the  ozone  photochemistry  schemes. 
In  the  lower  stratosphere  (100  to  lOhPa),  away  from  the 
ozone  hole,  biases  are  unconnected  with  the  photochemistry 
schemes.  Geer  et  al.  (2006a)  show  that  positive  biases  of 
up  to  20%  against  sonde  and  HALOE  are  likely  explained 
both  by  the  small  (~10%)  positive  bias  in  MIPAS  in  these 
regions,  and  by  poor  quality  transport,  a  known  deficiency 
in  stratospheric  data  assimilation  systems.  These  biases  are 
a  particular  problem  at  the  tropical  tropopause,  where  the 
DARC/Met  Office  analyses  are  50%  higher  than  ozoneson¬ 
des.  Similar  biases  were  found  in  many  ozone  analysis  sys¬ 
tems,  though  the  DARC/Met  Office  analyses  have  an  atyp- 
ically  large  bias  at  the  tropopause.  Biases  against  HALOE 
and  MIPAS  at  lOOhPa  and  below  should  be  treated  with 
much  caution  due  to  possible  cloud  contamination  of  the  ob¬ 
servations;  ozonesondes  are  much  more  reliable  here. 

We  also  examined  the  standard  deviations  of  difference  be¬ 
tween  analyses  and  observations.  The  results  were  in  general 
very  similar  to  those  seen  in  Geer  et  al.  (2006a),  and  typ¬ 
ical  of  many  ozone  data  assimilation  systems.  Significant 
differences  between  experiments  were  found  only  in  the  SH 
high  latitude  upper  stratosphere,  shown  in  Fig.  11.  HALOE 
and  MIPAS  have  different  sampling  patterns  (e.g.  Geer  et  al., 
2006a);  this  is  the  most  likely  reason  for  differences  between 
the  two  panels  of  Fig.  1 1 . 

The  following  sections  examine  these  biases  and  standard 
deviations  at  different  levels  in  the  atmosphere. 
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Fig.  8.  Mean  of  (analysis  -  HALOE)  ozone,  normalised  by  climatology,  in  latitude  bands  for  the  period  27  September  2003  to  5  November 
2003.  Vertical  scale  ranges  from  200  hPa  to  0.5  hPa.  See  colour  key  in  Fig.  3. 
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Fig.  9.  Mean  of  (analysis  -  ozonesonde)  ozone,  normalised  by  climatology,  in  latitude  bands  for  the  period  27  September  2003  to  5  November 
2003.  Vertical  scale  ranges  from  1000  hPa  to  10  hPa  (ozonesondes  do  not  usually  produce  data  above  this  level).  See  colour  key  in  Fig.  3. 
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Fig.  10.  Mean  of  (analysis  -  MIPAS)  ozone,  normalised  by  climatology,  in  latitude  bands  for  the  period  27  September  2003  to  5  November 
2003.  Vertical  scale  ranges  from  200  hPa  to  0.5  hPa.  See  colour  key  in  Fig.  3. 


4. 1  Upper  stratosphere  and  mesosphere 

In  the  upper  stratosphere  and  mesosphere,  the  LINOZ 
scheme  produces  negative  ozone  biases  that  reach  roughly 
40%  at  0.5  hPa  (Figs.  8  and  10).  This  is  a  known  prob¬ 
lem  with  the  scheme  (McCormack  et  al.,  2004;  Geer  et  al., 
2006a),  which  is  caused  by  an  excessive  net  loss  of  ozone 
driven  by  a  (P—L) o  term  that  is  up  to  several  orders  of 
magnitude  larger  than  that  in  the  other  schemes  in  this  re¬ 
gion  (Sect.  2.7,  Fig.  6).  Despite  such  large  differences, 
the  resulting  biases  are  no  more  than  ~40%  of  the  ozone 
field.  Why  should  this  be?  Equation  (7)  showed  that  the 
ozone  parametrization  causes  a  relaxation  to  steady  state  that 
is  particularly  strong  in  the  upper  stratosphere  and  meso¬ 
sphere.  In  Sect.  2.2  we  saw  that  an  error  e  in  the  ozone 
tendency  would  lead  to  a  steady  state  different  from  clima¬ 
tology:  Xs  v=X()+f r-  However,  because  the  photochemical 
relaxation  time  r  is  very  short  at  these  levels,  the  impact 
of  the  error  is  comparatively  small.  If  we  assume  the  other 
parametrizations  are  correct  then  we  can  estimate  the  error 
in  the  LINOZ  (P—L) o  term  from  Fig.  6  as  of  order  500% 
per  day  at  0.5  hPa.  From  Fig.  7,  the  corresponding  relaxation 
time  at  0.5  hPa  is  of  order  0.1  days.  Hence  we  would  expect 
an  error  in  ozone  of  around  50%,  in  agreement  with  Figs.  8 
and  10. 

There  are  positive  biases  in  Chem2D-OPP  v2.1  analyses 
in  the  upper  stratosphere  (1  hPa  to  lOhPa),  reaching  20%  in 


SH  high  latitudes,  compared  both  to  MIPAS  and  HALOE.  In 
the  same  region,  Chem2D-OPP  v2.1  shows  the  largest  stan¬ 
dard  deviations  against  MIPAS  and  HALOE  of  any  of  the 
analyses,  reaching  20%  against  MIPAS  (Fig.  11),  though  in 
other  latitude  bands  (not  shown)  there  is  little  difference  be¬ 
tween  experiments.  Figure  12  shows  examples  of  the  anal¬ 
ysed  ozone  fields  at  3.2  hPa.  It  appears  that,  at  latitudes  60°  S 
to  90°  S,  Cariolle  v2.1  represents  the  observed  ozone  field 
quite  realistically,  with  standard  deviations  of  ~6%  against 
HALOE  and  MIPAS  (Fig.  1 1).  In  contrast,  Cariolle  vl.O  and 
Chem2D-OPP  vO.  1  have  standard  deviations  of  ~10%,  sug¬ 
gesting  that  the  smaller  range  of  ozone  values  in  these  anal¬ 
yses  over  the  poles  (Fig.  12)  is  less  in  agreement  with  inde¬ 
pendent  data.  Structure  seen  in  the  ozone  field  over  the  pole 
in  the  Chem2D  v2.1  analyses  is  likely  erroneous,  given  the 
much  larger  (~20%)  standard  deviations. 

The  primary  explanation  for  the  erroneous  structure  in  the 
Chem2D-OPP  v2.1  ozone  fields  is  the  excessively  long  re¬ 
laxation  times  near  the  South  Pole  in  October  (shown  as  very 
low  sensitivities  in  Fig.  5h).  This  allows  the  temperature  term 
to  dominate,  but  the  temperature  term  in  Chem2D-OPP  v2.1 
also  has  problems.  We  have  already  seen  that  the  CIRA86 
temperature  climatology  used  with  Chem2D  is  substantially 
different  from  the  modelled  temperatures  (Fig.  2).  During 
the  vortex  breakup,  temperature  structures  are  often  far  from 
zonal  mean,  as  can  be  seen  from  Fig.  13.  Examination  of 
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HALOE  (201  obs) 


MIPAS  (5554  obs) 


Fig.  11.  Standard  deviation  of  (a)  (analysis  -  HALOE)  and  (b) 
(analysis  -  MIPAS)  ozone,  normalised  by  climatology,  in  the  lat¬ 
itude  band  90°  S  to  60°  S  for  the  period  27  September  2003  to 
5  November  2003.  See  colour  key  in  Fig.  3. 


the  individual  terms  shows  that  the  strongly  non-zonal  tem¬ 
perature  field  causes  the  temperature  term  in  Eq.  (1)  to  drive 
ozone  amounts  away  from  the  zonal  mean.  The  influence  of 
the  radiation  and  (P—L) o  terms  is  relatively  weak.  Oppo¬ 
sition  to  the  temperature  term  comes  mainly  from  the  ozone 
term,  which  returns  ozone  amounts  to  zonal  mean,  but  this  is 
erroneously  weak  (Fig.  5h).  Additionally,  these  experiments 
show  that  data  assimilation  increments  are  not  able  to  im¬ 
prove  the  ozone  field  where  the  temperature  term  is  so  strong: 
it  can  produce  changes  of  ozone  of  10%  in  a  day  with  just  a 
5K  perturbation  from  climatology  (see  Fig.  6c  above  5  hPa). 

McCormack  et  al.  (2006)  have  also  identified  prob¬ 
lems  with  the  temperature  term  in  model-only  runs  using 
Chem2D-OPP.  Particularly  in  the  polar  night,  where  the 
ozone  term  is  weak,  discrepancies  between  modelled  and 
CIRA86  climatology  temperatures  were  seen  to  cause  prob¬ 
lems. 

Cariolle  v2.1  shows  positive  biases  (Figs.  8  and  10)  at 
1  hPa  and  above,  reaching  a  maximum  of  20%  in  the  tropics 
at  0.5  hPa.  Biases  against  both  HALOE  and  MIPAS  are  typ¬ 
ically  smaller,  though  not  eliminated,  when  the  Fortuin  and 
Kelder  (1998)  climatology  is  used  instead  of  the  supplied  cli¬ 
matology.  This  suggests  that  the  climatology  supplied  with 
the  Cariolle  v2.1  scheme  is  slightly  biased  at  these  levels, 
and  that  replacing  it  with  the  Fortuin  and  Kelder  (1998)  cli¬ 
matology  can  remove  part  of  this  bias.  However,  even  after 
doing  this,  there  are  still  biases  of  order  10%  in  the  analyses 
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Fig.  12.  Ozone  fields  at  3.2  hPa  on  1  October  2003  from  (a) 
Cariolle  vl.0,  (b)  Cariolle  v2.1,  (c)  Chem2D-OPP  vO.l  and  (d) 
Chem2D-OPP  v2.1  analyses. 


versus  HALOE.  Section  2.7  has  shown  that  the  temperature 
and  ozone  terms  in  Cariolle  v2. 1  are  excessively  weak  in  the 
upper  stratosphere  and  mesosphere.  However,  the  relaxation 
time  is  still  roughly  a  day  (Fig.  7)  at  these  levels,  and  so 
it  still  causes  a  rapid  relaxation  to  climatology.  Following 
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the  arguments  in  Sect.  2.2,  it  is  likely  the  biases  could  be 
further  reduced  by  fine-tuning  the  temperature  climatology 
(To)  and  making  the  ozone  climatology  (xo)  consistent  with 
climatological  ozone  amounts  from  MIPAS.  However,  MI- 
PAS  and  HALOE  are  biased  by  ~5%  with  respect  to  each 
other  at  these  levels  (Geer  et  al.,  2006a),  so  even  if  the  model 
was  consistent  with  MIPAS,  there  would  necessarily  be  a 
bias  compared  to  HALOE.  Moreover,  this  would  improve  the 
simulation  for  the  period  studied  without  guarantee  of  its  ap¬ 
plicability  for  other  seasons. 

In  summary,  in  the  upper  stratosphere  and  mesosphere, 
LINOZ  is  unsuitable  for  use,  and  there  are  regional  biases 
in  the  Cariolle  v2.1  and  Chem2D-OPP  v2.1  experiments 
(respectively,  the  tropical  mesosphere  and  the  upper  strato¬ 
spheric  winter  vortex).  The  Cariolle  v2.1  biases  would  likely 
be  reduced  by  further  attention  to  the  temperature  and  ozone 
climatologies,  Tq  and  X(+  as  well  as  by  use  of  the  improved 
v2.8  coefficients  (Sect.  2.7).  Chem2D-OPP  v2.1  has  erro¬ 
neously  long  relaxation  times  in  the  south  polar  upper  strato¬ 
sphere  in  October,  but  this  problem  has  been  removed  in 
v2.6.  Excluding  LINOZ,  and  Cariolle  v2.1  and  Chem2D 
v2.1  in  the  problem  regions,  analyses  show  biases  in  the 
range  —10%  to  10%. 

4.2  Lower  stratosphere 

At  SH  high  latitudes  at  the  levels  where  ozone  depletion  takes 
place  in  the  ozone  hole  (100  to  40hPa),  the  Cariolle  vl.O  ex¬ 
periment  shows  roughly  20%  too  much  ozone  compared  to 
HALOE  and  sonde.  As  explained  in  more  detail  in  Geer  et  al. 
(2006b),  the  strong  radiation  term  of  the  vl.O  parametrization 
(Fig.  6d)  creates  erroneously  large  amounts  of  ozone  in  the 
ozone  hole.  The  other  experiments  show  positive  biases  of  no 
more  than  10%,  confirming  that  they  work  well  in  conjunc¬ 
tion  with  the  cold  tracer  heterogeneous  chemistry  scheme. 
The  biases  against  MIPAS  are  smaller  still,  indicating  that 
the  analyses  have  drawn  close  to  the  assimilated  MIPAS  ob¬ 
servations  and  the  remaining  biases  against  independent  data 
reflect  the  ~10%  positive  bias  between  MIPAS  and  indepen¬ 
dent  data  in  these  regions  (Geer  et  al.,  2006a). 

At  the  tropical  and  midlatitude  tropopause,  comparisons 
against  independent  data  show  no  difference  between  experi¬ 
ments.  This  is  despite  an  order  of  magnitude  difference  in  the 
sensitivity  of  the  ozone  terms  in  Cariolle  v2.1  and  Chem2D- 
OPP  v2.1  (Fig.  6b).  Photochemical  relaxation  times  are  in 
both  cases  extremely  long  (~100  days  and  ~1000  days,  re¬ 
spectively;  see  Fig.  7).  In  a  data  assimilation  system,  the 
chemistry  scheme  will  be  only  a  minor  part  of  the  ozone  bud¬ 
get  here,  which  will  instead  be  dominated  by  observational 
increments.  Hence,  evaluation  within  a  data  assimilation  sys¬ 
tem  is  not  able  to  distinguish  between  the  parametrizations 
or  to  suggest  which  may  be  more  correct.  Differences  would 
only  appear  in  relatively  long  free-model  runs,  which  would 
also  require  a  very  good  representation  of  tracer  transport. 


Fig.  13.  Analysed  temperature  (in  K)  at  3.2  hPa  on  1  October  2003. 

4.3  Troposphere 

No  ozone  observations  are  assimilated  below  approximately 
400  hPa,  and  the  assimilation  system  must  rely  instead  on 
modelled  ozone  photochemistry  and  transport  here.  We  have 
already  noted  that  the  photochemistry  parametrizations  are 
not  really  intended  for  use  in  the  troposphere.  However, 
many  current  ozone  assimilation  systems  do  use  Cariolle- 
type  parametrizations  in  the  troposphere  (Geer  et  al.,  2006a). 
It  is  worth  examining  how  large  are  the  resulting  errors. 

Figure  9  shows  biases  of  —50%  to  +100%  against  sonde 
in  the  troposphere.  The  largest  biases  tend  to  be  associated 
with  the  Cariolle  scheme:  Cariolle  vl.O  has  a  positive  bias  of 
over  100%  in  the  SH  polar  troposphere,  while  Cariolle  v2.1 
has  a  positive  bias  of  80%  in  the  tropical  lower  troposphere. 
The  latter  bias  has  also  been  seen  in  assimilation  runs  with 
the  MOCAGE/PALM  data  assimilation  system  using  Cari¬ 
olle  v2.1  photochemistry  (Geer  et  al.,  2006a).  By  a  detailed 
examination  of  the  coefficients  (not  shown),  it  appears  that 
these  positive  biases  arise  because  of  ozone  production  by 
the  (P—L) o  term,  which  is  substantially  larger  than  those  in 
LINOZ  or  Chem2D-OPP  in  the  lower  troposphere.  Figure  6a 
gives  some  indication  of  this.  This  model  bias  in  the  equa¬ 
torial  lower  troposphere  has  been  corrected  from  v2.3  of  the 
Cariolle  scheme  onwards. 

Other  schemes  are,  because  of  the  way  they  have  been  im¬ 
plemented  in  these  experiments,  dominated  by  the  relaxation 
to  climatology  of  the  ozone  term.  For  these  schemes,  biases 
are  in  general  within  —20%  to  +20%,  except  in  the  NH  polar 
lower  troposphere,  where  biases  can  be  as  large  as  —40%, 
though  results  are  based  on  very  few  sondes.  The  correct 
handling  of  tropospheric  ozone  remains  an  outstanding  chal¬ 
lenge  in  ozone  data  assimilation  and  a  fast  parametrization 
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of  tropospheric  ozone  chemistry  is  not  yet  available.  Hence, 
a  pragmatic  temporary  solution  is  to  relax  ozone  amounts  to 
climatology.  In  these  experiments  it  is  a  2-D  climatology,  but 
use  of  a  3-D  one  like  Logan  (1999)  would  capture  the  main 
zonal  features  of  the  tropospheric  ozone  distribution,  which 
are  visible  in  total  column  observations  in  the  tropics.  This 
would  be  important  if  total  column  observations  were  assim¬ 
ilated,  so  as  to  avoid  the  aliasing  of  tropospheric  biases  into 
changes  in  stratospheric  ozone. 


5  Conclusion 

This  study  has  examined  ozone  analyses  from  the  Met  Office 
stratosphere/troposphere  data  assimilation  system.  MIPAS 
ozone  retrievals  were  assimilated  for  the  period  23  Septem¬ 
ber  2003  until  5  November  2003  into  a  set  of  experi¬ 
ments,  each  using  a  different  linear  ozone  photochemistry 
parametrization.  Heterogeneous  chemistry  was  parametrized 
using  a  cold  tracer  scheme  (Eskes  et  ah,  2003),  and  remained 
fixed  throughout  the  experiments.  Analysed  ozone  was  val¬ 
idated  against  independent  observations  from  HALOE  and 
ozonesondes.  None  of  the  ozone  parametrizations  is  specifi¬ 
cally  intended  for  use  in  the  troposphere,  where  modelling 
and  observing  ozone  remains  a  substantial  challenge,  but 
results  were  in  general  good  in  the  stratosphere  and  lower 
mesosphere. 

We  have  also  introduced  a  simple  offline  method  to  com¬ 
pare  the  effect  of  the  coefficients  in  different  parametriza¬ 
tions.  This  is  done  by  calculating  the  instantaneous  rate  of 
change  of  ozone  that  would  result  from  a  representative  per¬ 
turbation  in  ozone  or  temperature.  With  this  technique  we 
were  able  to  learn  a  lot  about  the  different  schemes  and  to 
show  that  our  results  can  be  generalised  beyond  the  relatively 
short  period  of  the  assimilation  experiments. 

In  the  upper  stratosphere  and  mesosphere,  outside  the  po¬ 
lar  night,  the  photochemical  lifetime  of  ozone  is  short,  and 
analyses  are  controlled  mostly  by  the  ozone  parametrization, 
not  by  the  assimilated  observations.  Here,  the  parametriza¬ 
tion  causes  modelled  ozone  to  relax  quickly  to  a  steady  state. 
Hence,  even  though  rate  coefficients  may  vary  between  the 
different  schemes,  often  the  most  important  factor  in  con¬ 
trolling  the  modelled  or  assimilated  ozone  distribution  is  ac¬ 
tually  the  temperature  and  ozone  climatology  supplied  with 
each  scheme.  This  paper  illustrates  the  problems  that  can 
be  caused  when  there  are  discrepancies  between  the  clima¬ 
tologies  of  temperature  and  ozone  used  within  the  ozone 
parametrization  and  their  equivalents  in  the  GCM  and  in  the 
assimilated  observations.  Improvements  could  be  made  by 
making  the  climatologies  of  ozone  and  temperature  in  the 
schemes  more  consistent  with  their  equivalents  in  the  GCM. 
In  data  assimilation  applications,  models  should  also  be  unbi¬ 
ased  with  respect  to  observations.  In  regions  where  the  pho¬ 
tochemical  lifetime  is  short,  this  means  the  schemes  should  if 
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possible  use  an  ozone  climatology  based  on  the  assimilated 
data. 

In  the  mid  and  lower  stratosphere,  analyses  strongly  ben- 
efitted  from  the  assimilation  of  good  quality  MIPAS  data 
with  a  relatively  high  vertical  resolution.  Photochemical 
timescales  are  relatively  slow  here  (>100  days);  any  slow- 
growing  model  biases  would  be  corrected  by  the  observation 
increments,  which  are  added  on  much  shorter  timescales. 
However,  a  system  with  more  sparse  observational  coverage 
for  ozone  (typical  of  operational  systems)  might  reveal  larger 
biases  and  differences  between  the  photochemistry  schemes. 
It  is  clear,  however,  that  the  data  assimilation  method  is  less 
able  to  test  the  photochemistry  parametrizations  in  regions 
where  their  influence  is  weaker  such  as  the  lower  strato¬ 
sphere. 

A  particular  example  is  the  differing  relaxation  times  (r) 
in  the  Cariolle  v2.1  and  Chem2D-OPP  v2.1  schemes  in  the 
lower  stratosphere.  Here,  Chem2D-OPP  v2.1  is  roughly  an 
order  of  magnitude  less  sensitive  to  ozone  variations  than 
Cariolle  v2.1,  mainly  because  Chem2D-OPP  v2.1  neglects 
Brx  chemistry,  but  also  likely  because  of  the  different  cal¬ 
culation  methods  used.  However,  despite  the  order  of  mag¬ 
nitude  differences  between  the  parametrizations,  we  saw  no 
difference  in  analysed  ozone  in  the  lower  stratosphere. 

We  see  that  data  assimilation  is  best  for  testing  the  photo¬ 
chemistry  schemes  in  regions  where  they  have  a  more  dom¬ 
inant  control  on  the  ozone  amounts,  such  as  in  the  mid  and 
upper-stratosphere.  Data  assimilation  does  however  force  the 
schemes  to  operate  in  an  environment  as  close  to  reality  as 
can  be  provided,  which  model  free-runs  cannot  do.  How¬ 
ever,  it  is  clear  that  the  results  in  this  paper  would  be  use¬ 
fully  complemented  by  long-term  GCM  or  CTM  runs  which 
would  allow  slow-growing  problems  to  be  identified,  and  by 
comparisons  to  more  detailed  chemistry  schemes.  Chem2D- 
OPP  v2  and  Cariolle  v2  have  been  separately  tested  in  long- 
period  model  runs  by  McCormack  et  al.  (2006)  and  Cariolle 
and  Teyssedre  (2007). 

We  next  summarise  the  results  for  each  parametrization 
in  turn.  At  the  time  the  assimilation  runs  were  made,  we 
only  had  available  the  LINOZ,  Cariolle  vl.O  and  v2.1,  and 
CHEM2D-OPP  vO.l  and  v2.1  parametrizations.  Before  final 
publication,  we  have  been  able  to  examine  the  latest  available 
parametrizations  Cariolle  v2.8  and  CHEM2D-OPP  v2.6,  us¬ 
ing  the  offline  method  of  Sect.  2.7.  A  number  of  the  dis¬ 
crepancies  we  identified  in  Cariolle  v2.1  and  CHEM2D-OPP 
v2. 1  have  now  been  fixed,  partly  in  response  to  this  study.  In 
this  work  we  have  also  examined  a  number  of  problems  (such 
as  that  of  LINOZ  in  the  upper  stratosphere)  that  are  already 
known,  but  which  are  presented  here  in  a  wider  context.  In 
the  following  summary,  reference  is  made  to  the  earlier  work 
in  such  cases.  We  found: 

-  LINOZ  is  unsuitable  for  use  above  lOhPa  (McCormack 
et  al.,  2004;  Geer  et  al.,  2006a),  though  below  this  level, 
it  generally  works  as  well  as  the  other  schemes. 
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-  Cariolle  vl.O  has  an  excessively  strong  dependence  on 
overlying  column  ozone  in  the  radiation  term  (Geer 
et  al.,  2006b),  which  causes  excessive  ozone  production 
in  the  ozone  hole,  leading  to  ozone  values  20%  higher 
than  sonde  in  the  SH  lower  stratosphere  (100  to  30hPa) 
in  the  current  experiments.  There  are  also  unrealisti¬ 
cally  high  ozone  values  near  the  surface  around  Antarc¬ 
tica.  Since  these  coefficients  were  generated  many  years 
ago  (Cariolle  and  Deque,  1986),  it  is  hard  to  identify  ex¬ 
actly  where  these  problems  may  have  come  from.  Else¬ 
where,  Cariolle  vl.O  works  in  general  as  well  as  the 
other  schemes. 

-  The  Cariolle  v2.1  scheme  performed  well  in  general.  A 
~20%  positive  bias  in  the  tropical  upper  stratosphere 
and  mesosphere  was  caused  largely  by  the  supplied 
ozone  climatology,  and  could  be  partially  corrected  by 
substituting  the  Fortuin  and  Kelder  (1998)  climatology. 
Temperature  and  ozone  terms  were  also  too  weak  in  the 
upper  stratosphere  and  mesosphere,  and  this  has  been 
corrected  in  later  versions.  There  was  also  an  80% 
positive  bias  in  the  tropical  lower  troposphere,  linked 
to  excessive  ozone  production  in  the  (P—L) o  term, 
which  has  been  largely  corrected  in  later  versions  of  the 
scheme.  Cariolle  v2.1  analyses  had  particular  success  in 
the  SH  high  latitude  upper  stratosphere  (6  hPa  to  1  hPa), 
which  experienced  strong  minor  warmings  through  Oc¬ 
tober  2003,  with  synoptic  temperature  variations  of 
~50K.  Here,  Cariolle  v2.1  analyses  showed  only  6% 
standard  deviation  of  difference  against  HALOE,  com¬ 
pared  to  10%  for  most  other  schemes  and  20%  for 
Chem2D-OPP  v2.1. 

-  Analyses  using  Chem2D-OPP  vO.l  (only  the  first  two 
terms  in  Eq.  1)  and  v2.1  (all  terms)  performed  well  in 
general.  Chem2D-OPP  v2.1  produced  unrealistic  ozone 
structures  in  the  SH  polar  upper  stratosphere  (6  hPa  to 
1  hPa),  due  primarily  to  an  error  in  the  ozone  relaxation 
times,  which  has  been  fixed  from  v2.4  onwards.  The 
problem  was  also  likely  exacerbated  by  discrepancies 
between  analysed  temperatures  and  the  CIRA86  tem¬ 
perature  climatology  supplied  with  the  scheme.  A  new 
temperature  climatology  is  currently  being  tested,  and 
shows  promising  results. 

These  tests  are  based  on  a  short  period  in  September  to 
November.  The  precise  behaviour,  for  example  in  terms  of 
the  bias  in  any  scheme,  will  likely  vary  from  region  to  re¬ 
gion  and  through  the  year.  Also,  the  exact  sizes  of  the  bi¬ 
ases  are  specific  to  our  data  assimilation  system  and  the  MI- 
PAS  observations.  Nevertheless,  we  believe  that  these  re¬ 
sults  have  wider  significance.  First,  most  of  the  problems  we 
have  illustrated  can  be  explained  by  specific  features  of  the 
coefficients  in  the  various  schemes.  These  features  are  typi¬ 
cally  persistent  throughout  the  year  and  at  different  latitudes 
(Figs.  4  and  5).  Hence,  these  figures  can  be  used  to  predict 


where  similar  behaviour  is  likely  to  occur  at  other  times  of 
the  year.  Second,  where  the  influence  of  the  photochemistry 
schemes  is  strong,  i.e.  in  the  mid  and  upper-stratosphere,  all 
models  and  data  assimilation  systems  are  likely  to  be  affected 
in  ways  similar  to  those  seen  here.  For  example,  in  experi¬ 
ments  using  the  KNMI  system  with  the  assimilation  of  only 
total  column  ozone  observations  (Geer  et  al.,  2006a),  the  use 
of  LINOZ  caused  similar  biases  in  the  upper  stratosphere  to 
those  seen  in  Fig.  10. 

Future  developments  of  the  ozone  photochemistry 
schemes  should  include  a  better  treatment  of  the  troposphere. 
Also,  none  of  these  schemes  simulates  the  diurnal  cycle  of 
ozone  above  0.5  hPa,  yet  this  has  a  strong  influence  on  mod¬ 
elled  temperatures  (Sassi  et  al.,  2005). 
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